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System-size dependence of strangeness production 
in nucleus-nucleus collisions at ^/snn = 17.3 GeV 
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Emission of tt , K , (f> and A was measured in near-central C+C and Si+Si collisions at 158 AGeV 
beam energy. Together with earlier data for p+p, S+S and Pb+Pb, the system-size dependence 
of relative strangeness production in nucleus-nucleus collisions is obtained. Its fast rise and the 
saturation observed at about 60 participating nucleons can be understood as onset of the formation 
of coherent systems of increasing size. 

PACS numbers: 25.75.-q 



Most recent experimental studies of particle produc- 
tion in relativistic nucleus-nucleus collisions have fo- 
cussed on heavy reaction partners because of the expecta- 
tion that the transition to a deconfined state of strongly 
interacting matter [ij is more likely to occur in large col- 
lision systems. However, earlier experiments with lighter 
beams at the CERN SPS accelerator had demonstrated 
that already in '^^S-l-S collisions |2] relative strangeness 
production — which has been discussed as a possible indi- 
cator for such a transition 'sf — is significantly enhanced 



compared to p+p collisions and, as found later, is close 
to the situation in ^ospb-hPb |10. 

In order to understand the origin of this enhancement 
a systematic study of its dependence on the size of the 
collision system is required. In this Letter we discuss 
systems with similar geometry, i. e. near-central and ap- 
proximately symmetric A-|- A collisions between nuclei of 
varying mass number A. New data on the production of 
TT^, K^, (/>- mesons and A-baryons will be presented for 
the systems C-l-C and Si+Si. They were obtained with 



the NA49 spectrometer at the CERN SPS at 158 AGeV 
beam energy, corresponding to a center-of-mass energy 
of y/sNN = 17.3 GeV. The results will be discussed to- 
gether with available data for p+p and Pb+Pb at 158 
AGeV and S+S at 200 AGeV beam energy. It will be 
shown that strangeness enhancement develops fast with 
increasing size of the collision system reaching the level of 
Pb+Pb already at about 60 participating nucleons. We 
argue that this result can be understood assuming the 
formation of coherent systems of increasing size. 

NA49 is a fixed-target experiment at CERN |a| using 
external SPS beams. For the present study a fragmented 
Pb beam of 158 AGeV was used |y]. The fragments were 
selected by magnetic rigidity {Z/A = 0.5) and by spe- 
cific energy loss in transmission detectors in the beam 
line. The "C-beam" as defined by the online and offline 
trigger was a mixture of ions with Z=6 and 7 (inten- 
sity ratio 69:31) in the case of the meson measurements, 
and pure Z=6 ions for the A. As "Si-beam" we used ions 
with Z=13,14, and 15 (35:41:24). 3 mm and 10 mm thick 
carbon targets (2.4% and 7.9% interaction length) were 
chosen, and a 5 mm thick Si target (4.4%). The veto 
calorimeter placed 20 m downstream of the target ac- 
cepts nearly all beam particles, projectile fragments and 
spectator neutrons and protons j^ . By setting an upper 
threshold for the energy deposit, the (15.3 ± 2.4)% most 
central C-l-C and (12.2 ± 1.8)% most central Si-|-Si in- 
teractions were selected [3. After all cuts typically 45k 
events were analyzed for the mesons and 200k for the 
A. Simulations based on the VENUS model of nucleus- 
nucleus interactions 8] allow to calculate impact param- 
eter range and mean number of wounded nucleons '9|: 
b < 2.0 fm for C-hC and b < 2.6 fm for Si+Si, for 
(A'wound) see table |n] In A -I- A reactions some nucleons 
in the collision may participate only through secondary 
cascading processes. In this paper, we do not include 
such nucleons in counting the number of wounded nucle- 
ons iVwound, but we do include them in the number of 
participating nucleons iVpart, which is also provided by 
VENUS. These definitions differ from those assumed for 
Glauber calculations. For the Apart values given in ta- 
ble |nl the numbers from VENUS are averaged with those 
derived from the measured energy of the projectile spec- 
tators in the veto calorimeter. 

Charged particles produced in beam-target inter- 
actions were measured in the NA49 large-acceptance 
hadron spectrometer consisting of magnets and time pro- 
jection chambers (TPCs), two of them immediately be- 
hind the target and inside the magnetic field, two addi- 
tional ones further downstream on both sides of the beam 
(for details see Q). By tracking and measurement of the 
specific energy loss {dE/dx) in the TPCs the momen- 
tum and the identity of the particles were determined. 
In the (lab)momentum range 4 < p < 50 GeV/c and 
with a relative dE/dx resolution ranging from about (3- 
6)% depending on track length and momentum 6] pi- 



ous, kaons and protons could be resolved on a statisti- 
cal basis by decomposing the dE/dx spectra in (j>,pt) 
bins which are then transformed to {y,pt) or (y, mt) bins 
(pt — transverse momentum, mt = transverse mass, y 
= cm. rapidity). Within the geometrical acceptance of 
the spectrometer which covers most of the forward hemi- 
sphere in the center-of-mass system of the collision pious 
and kaons were selected in a range of the azimuthal an- 
gle (about 28% of 2tt) where small background and a flat 
acceptance is ensured. Acceptance and corrections for in- 
flight decay of kaons (about 2% losses in addition) were 
determined by GEANT [lOj simulations. The contribu- 
tion of pions from weak decays of hadrons (mainly K|g) is 
subtracted using events from the VENUS model (about 
5%). Care is taken that the VENUS multiplicities agree 
with data, otherwise they are scaled accordingly. The (p- 
niesons were measured via their strong decay into K+K^ 
pairs (branching ratio 49.1%). Yields were derived in 
y- or pt-bins from the corresponding peak in the K+K^ 
invariant-mass spectra after subtracting the combinato- 
rial background; the overall acceptance in the forward 
hemisphere amounts to about 45%. The technique is de- 
scribed in jllj ; for more results also on the other mesons 
see |2|. A-hyperons decay weakly via A — > pvr^ (branch- 
ing ratio 63.9%). They were identified on the basis of 
their typical VO decay topology and the p7r~ invariant 
mass; for the method see |5[. The P(-integrated com- 
bination of acceptance and detection efficiency ranges 
from 8% in the most backward y-hin to about 57% at 
midrapidity. Since all S*' decay electromagnetically to 
A-baryons, A denotes the sum of both. For Si+Si an addi- 
tional correction of (5±3)% for the reduced A reconstruc- 
tion efficiency was applied due to the higher track den- 
sity in these reactions. GEANT simulations showed that, 
varying slightly with pt and y, 90% of the A-baryons from 
weak decays of S~ and S" contribute to the measured A- 
yield. Using yield estimates from statistical model fits 
|l2l | a feeddown correction of (9±3)% and (10±3)% was 
applied for A from C+C and Si+Si, respectively. 

For all particles under study yields were obtained in 
{y,Pt) bins. The transverse distributions, those close to 
midrapidity being shown in fig.Q were fitted with a ther- 
mal ansatz d^ N / (rritdmtdy) ex exp(— mt/T). The fit was 
used to extrapolate to pt-regions not covered by the mea- 
surement, and the yield dN/dy was obtained by summing 
the measured bins and this extrapolation. The latter cor- 
rection was of the order of a few percent typically. An 
exception were pions at midrapidity (—0.13 < y < 0.33) 
for which the correction factors became as large as 2 to 
3 due to the limit p < 4 GeV/c for particle identifica- 
tion. In this case, spectral shapes from the NA35 S+S 
[l3j measurements were used, this way taking the known 
low-pt pion enhancement into account. In order to es- 
timate the systematic error in this phase space region 
different analysis strategies were investigated: When ex- 
trapolating the pion yields with a straight exponential 
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FIG. 1; Transverse mass distributions for n and K in 0.1 < 
y < 0.33, for (p averaged over < y < 1.8, and for A in 
\y\ < 0.4 measured in C+C (•) and Si+Si (a) collisions. Only 
statistical errors are shown. 



with an inverse slope T as given in table Q] yields would 
be lower by (6-9)%. An h~ analysis as performed in |j| 
with corrections for K^ and p gives (5-9)% higher yields. 
This uncertainty mainly determines the systematic er- 
ror on the midrapidity and total pion yields. The pion 
rapidity distributions dN/dy (fig. ^ were approximated 
by a superposition of two Gaussians displaced symmet- 
rically around mid-rapidity by a shift iyo- For kaons 
and the (/)-meson a single Gaussian was used. Integration 
gives the total average yield (TV) of each particle. To 
obtain the yields of A-baryons their y-distribution was 
fitted by a shape adopted from p+p and S+S measure- 
ments. The mean of both approximations was taken for 
C+C; for Si+Si only the shape from S+S was used (lines 
in fig. 12)1. For the (/)- meson with its much lower abun- 
dance the pt-distribution could only be obtained when 
averaged over a wider rapidity range, < y < 1.8. The 
values for dN/dy were determined in j/-slices using the pt- 
range < pt < 1.5 GeV/c and the averaged T-parameter 
for extrapolation (fig. O. 

Parameters for the kinetic distributions of all studied 
particles are given in tabled their yields in table Hll 

The strange-particle yields were normalized to the 
mean number of charged pions per event, (tt*) = ((7r+) + 
(7r~))/2. In this way they represent the relative contribu- 
tion of strange quarks in the final state. The results are 
shown in fig.|31together with data for minimum-bias p+p 
[Ellllllll and Pb+Pb interactions (5% most central, A- 
baryons 10% most central) |j, y, llll taken at the same 
energy. The A-yield is scaled by the ratio of the num- 
ber of wounded nucleons to account for the difference in 
the centrality selection, and the correction for feeddown 
from S-baryons is applied pj. Also plotted are the data 
for S+S (2% most central) obtained at the slightly higher 
energy of ^s^vw = 19.4 GeV [i|. A-yields for S+S were 
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FIG. 2: Rapidity distributions for C+C (•) and Si+Si (A) 
collisions. Open symbols are reflected at midrapidity. Error 
bars are statistical, systematic errors at midrapidity are indi- 
cated by the rectangular boxes. The curves represent fits to 
the data as described in the text. 
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FIG. 3: Experimental ratios of {K+), (K"), (<^), and (A) to 
(tt ) plotted as a function of system size (t p+p, • C+C and 
Si+Si, A S+S, ■ Pb+Pb). Statistical errors are shown as error 
bars, systematic errors if available as rectangular boxes. The 
curves are shown to guide the eye and represent a functional 
form a- b- exp(-{7Vpart)/40). At (TVpart) = 60 they rise to 
about 80% of the difference of the ratios between A'^part = 2 
and 400. 



c+c 


T [MeV] 


dN/dy 


ay 


yo 


t:+ 


171 ± 10 


5.6 ± 0.6 


1.12 ± 0.09 


0.87 ± 0.07 


TT~ 


171 ± 10 


5.7 ± 0.6 


1.11 ± 0.09 


0.85 ± 0.07 


K+ 


188 ± 10 


0.85 ± 0.085 


1.23 ± 0.08 




K- 


185 ± 10 


0.53 ± 0.053 


1.10 ± 0.09 




(p 


189 ± 28 


0.062 ± 0.008 


1.16 ± 0.10 




A 


199 ± 15 


0.24 ± 0.038 






Si+Si 


r [MeV] 


dN/dy 


ay 


yo 


7T + 


173 ± 10 


14.8 ± 1.5 


1.04 ± 0.08 


0.91 ± 0.07 


TT~ 


178 ± 10 


15.0 ± 1.5 


1.05 ± 0.08 


0.89 ± 0.07 


K+ 


192 ± 10 


2.4 ± 0.24 


1.24 ± 0.08 




K" 


196 ± 10 


1.5 ± 0.15 


1.15 ± 0.07 







220 ± 28 


0.19 ± 0.02 


1.27 ± 0.10 




A 


235 ± 16 


0.88 ± 0.13 







TABLE I: Parameters of the measured kinetic distributions. 
For pions and kaons T holds for 0.1 < y < 0.33, dN/dy for 
—0.13 < y < 0.33; for A both values are given for \y\ < 0.4. 
For the (/f>-meson T is determined in < j/ < 1.8, dN/dy in < 
y < 0.4. Only the dominating errors are given, i.e. statistical 
errors for the and systematic errors for the other hadrons. 
Statistical errors of the latter are less than half the systematic 
error; for the (p the systematic error is either the same (dN/dy) 
or half of the statistical error. 



C+C 



Si+Si 



\^^-wound/ 


14 + 2 


37 + 3 


(A^part) 


16.3 + 1 


41.4 + 2 


(-+) 


22.4 + 0.3 + 1.6 


56.6 + 0.7 + 4 


^-) 


22.2 + 0.3 + 1.6 


57.6 + 0.6 + 4 


(K+) 


2.54 ± 0.03 ± 0.25 


7.44 + 0.08 + 0.74 


(K-) 


1.49 ± 0.05 ± 0.15 


4.42 + 0.04 + 0.44 


(0) 


0.18 + 0.01 + 0.02 


0.66 + 0.03 + 0.08 


(A) 


1.32 + 0.05 + 0.32 


3.88 + 0.16 + 0.56 



TABLE II: Particle yields for C+C (15.3% most central col- 
lisions) and Si+Si (12.2%) at y/s^ = 17.3 GeV. The first 
error is statistical, the second systematic. 



scaled down by 10% to account for feeddown from S- 
baryons assuming the same S/A ratio as in Si+Si and 
using the information on the magnitude of the feeddown 
contribution given in J2| . 

Using the mean number of nucleons participating in 
the collision ((A'^part)) as a measure of system size the 
resulting dependence of the relative strangeness pro- 
duction (fig. |21l is characterized by a fast rise at small 
(.^part) values, reaching the level of Pb+Pb interactions 
at (A'part) ~ 60. Note that the same shape of the 
enhancement curve holds for all measured strangeness- 
carrying particles including the (/)-meson with only hid- 
den strangeness. 

The system-size dependence of relative strangeness 
production indicates basic changes in the reaction mech- 



anism due to the transition from an isolated N+N inter- 
action to the situation in A-\- A where nucleons undergo 
several subsequent collisions within short time intervals 
and where many of these collision sequences occur close 
in space. Microscopic transport models such as UrQMD 
|l5j calculate position and time of each collision in the 
center-of-mass system of the reaction. This allows to ex- 
tract a collision density per unit space and time interval 
averaged over the interpenetration phase of the colliding 
nuclei. This collision density is found to increase with 
the size of the colliding nuclei; for central Si+Si reac- 
tions about 4-5 interactions per fm'^ and fm/c are reached 
[3, |l6j . It is reasonable to assume that neither the sub- 
sequent collisions nor the decay of the created objects 
proceed independently. We propose that the reaction 
volume should thus become quantum-mechanically co- 
herent. Apparently, these conditions influence the global 
strangeness production rate. 

Statistical models for hadron production (e.g. |l7lll8l|) 
are very successful in reproducing the particle ratios in 
p+p as well as in A + A collisions. In these models 
the implementation of strangeness conservation results 
in a strong volume dependence of strangeness produc- 
tion |19ll2rll2ll |. It is natural to invoke this effect, known 
as canonical strangeness suppression, as a qualitative ex- 
planation for the data shown in fig. |31 Quantitatively, 
however, discrepancies remain. In the model of Tounsi 
and Redlich |21|, based on a hadronic scenario, as well as 
in the model of Rafelski and Danos [lS| , based on both, a 
partonic and a hadronic scenario, the strangeness satura- 
tion point (as defined in the caption of fig. 01 for jS*! = 1 
particles is close to iVpart — 10, in contrast to the exper- 
imental value of A^part = 60. One can reconcile theory 
with experiment assuming that only parts of the reaction 
volume form the statistical ensembles used in the mod- 
els. This implies a modification of the linear relationship 
between volume and number of participants which was 
assumed so far. 

To summarize, the data obtained for near-central sym- 
metric nucleus-nucleus collisions with increasing number 
of participating nucleons show a steep rise of the relative 
strangeness content in the produced hadrons reaching a 
saturation value at around (A'^part) '^ 60. It is argued 
that this behavior is correlated with the formation of ex- 
tended systems of dense and excited matter which collec- 
tively decay to the observed hadrons. Within the statisti- 
cal model for hadron production this would qualitatively 
explain the increase of relative strangeness production, 
because the constraints due to strangeness conservation 
are diminishing with increasing size of the hadronizing 
systems. Strangeness saturation thus occurs when the 
grand-canonical limit can be applied |l9ll2Cll2lLl22| . The 
experiment shows that for singly strange particles and the 
0-meson at y^JJTN = 17.3 GeV this happens already in 
central collisions between nuclei with A about 30. 
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